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A new formula is presented for the calculation of matrix elements between multi-quasiparticle 
Hartree-Fock-Bogoliubov (HFB) states. The formula is expressed in terms of the Pfaffian, and is 
derived by using the Fermion coherent states with Grassmann numbers. It turns out that the for- 
mula corresponds to an extension of generalized Wick's theorem and simplifies the combinatorial 
complexity resulting from practical applications of generalized Wick's theorem by unifying the tran- 
sition density and the transition pairing tensor in the HFB theory. The resultant formula is simpler 
and more compact than the traditional description of matrix elements of general many-body oper- 
ators. In addition, through the derivation of our new formula, we found that the Pfaffian version 
of the Lewis Carroll formula corresponds to the relation conjectrured by Balian and Brezin for the 
HFB theory in 1969. 
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I. INTRODUCTION 

A Fermion coherent state with Grassmann numbers 
has been useful not only to particle physics, but also 
to quantum many-body physics [![. For instance, tradi- 
tional perturbation theories both for zero and finite tem- 
perature quantum systems can be formulated by means 
of the functional integral with the Fermion coherent state 
and Grassmann numbers [l|. With this approach, we can 
obtain a new insight about Green's functions, partition 
functions in statistical mechanics, and many-body matrix 
elements (overlaps for many-body states), and so on. 

Such a Fermion coherent state with Grassmann num- 
bers was recently applied to nuclear structure physics by 
Robledo Q. His discovery is that the Thouless form of 
Hartree-Fock-Bogoliubov (HFB) state can be mapped to 
the Fermion coherent state with the help of Grassmann 
numbers. It sheds a new light on hitherto the established 
formulation and the long-standing unresolved issues in 
nuclear many-body systems. 

In nuclear structure physics, the HFB method has 
played a central role in describing interacting nucleons 
through a short-range effective force (e.g., pairing inter- 
action) as well as a long-range force (e.g., quadrupole 
interaction) Q. With various constrained conditions, 
the HFB has succeeded in giving a good mean-field de- 
scription for various many-body states. However, in an 
attempt to go beyond the mean-field description, espe- 
cially in the case of three-dimensional angular momen- 
tum projection, there has been a difficulty originating 
from the long-standing problem in the phase determi- 
nation of norm-overlap kernels through the Onishi for- 
mula |4J . This problem has been thoroughly investigated 
by many authors, for instance, in Refs. 4-8]. Many of 
them relied upon the analytic continuity approach for the 
phase determination, which can be carried out with the 
Onishi formula |4J. It was demonstrated recently that this 
approach works even for high-spin states in the presence 
of many nodal lines @ . 

On the contrary to these pursuits, Robledo was suc- 
cessful in deriving a new formula for the HFB norm- 



overlap kernel in terms of the Pfaffian. The formula con- 
tains no ambiguity in phase. This achievement was made 
by means of the Fermion coherent states with Grassmann 
numbers |2|. By his formula, the perennial problem has 
been finally solved in a mathematically elegant fashion. 
We extended his methodology to find another formula, 
which is app licable both to even- and odd-mass nucleon 
systems [Tol ]. 

A further extension should be to derive a formula ap- 
plied to multi-quasiparticle states. A standard approach 
to evaluate matrix elements between multi-quasiparticle 
states is to use generalized Wick's theorem. By this the- 
orem, evaluation of the matrix elements between multi- 
quasiparticle states is reduced to a combinatorial problem 
with contractions, which correspond to the density ma- 
trix and the pairing tensor in the HFB theory (the transi- 
tion density and the transition pairing tensor in the case 
of projected HFB). However, due to the combinatorial 
nature, the number of combinations of the contractions 
rapidly increases as we consider more complicated multi- 
quasiparticle states. Therefore, it is important to derive 
a systematic, compact and closed form for the formula in 
order to calculate such matrix elements effectively. Perez- 
Martin and Robledo derived a formula based on the finite 
temperature formalism [11]. Recently, different formulae 
expressed in terms of the Pfaffian have been propose d by 
Bertsch and Robledo [ill], and by Avez and Bender (l3j. 
although the latter does not correspond to generalized 
Wick's theorem. 

In this paper, we also present a new compact and 
closed form (but more compact than the above former 
two) with the Pfaffian by means of the Fermion coherent 
state with Grassmann numbers. In Ref . flfjl] . we have al- 
ready presented a Pfaffian formula to calculate the norm- 
overlap kernels of one-quasiparticle states. This formula 
is extended to multi-quasiparticle cases in the present 
work, by introducing a generating function expressed 
by Grassmann numbers. It turns out that the formula 
corresponds to an extension of generalized Wick's the- 
orem. By comparing our formula and the one by Avez 
and Bender [13(, we can give an affirmative answer to the 
Balian-Brezin conjecture [l4| . It is surprising to find that 
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a usefulness of the Pfaffian in the HFB theory has been 
already suggested in the their paper in 1969. 

The present paper is organized as follows. Section II 
is devoted to an explanation of the general mathemati- 
cal structure of matrix elements of many-body operators. 
Section III is for the Grassmann representation of the ma- 
trix elements. Section IV is a brief review on the Grass- 
mann integrations. Section V and VI present a derivation 
of a Pfaffian form of the considered matrix elements, and 
a discussion concerning generalized Wick's theorem. In 
Section VII, a discussion is given on the Balian-Brcsin 
conjecture. Section VIII is the summary of the present 
paper. 

Because the references necessary for the present paper 
cover a wide range of fields in phyiscs and mathematics, 
such as particle physics, mathematical physics, combina- 
torics in pure mathematics and so on, it is convenient to 
structure this paper in a self-contained manner. 



II. MATRIX ELEMENT FOR 
MULTI-QUASIPARTICLE STATES 

In this paper, the central representation of the HFB 
wave function is the Thouless formQ, which is presented 
by 




(1) 



where p takes or 1, and Af is a normalization factor. 
The creation (annihilation) operator of a single-particle 
orbit i is expressed as cj (ej). The total dimension of the 
single-particle model space is M. 0) is the vacuum state 
to this basis, that is, Cj|0) = 0. Anti-symmetric matrix 
Z is expressed in terms of the Bogoliubov transformation 
(or the so-called U and V matrices), as 

Z {p) = (7(p)(^(p))-1)*, (2) 
An m-quasiparticle HFB state is given by 

l^ 1 ) ,., fe J=# + ---/3£ + l^ ) ), (3) 

where fci, fe, • - • , k m are indices for quasiparticle states, 
and /3's (/3^'s) are quasiparticle annihilation (creation) 
operators, which are defined with the Bogoliubov trans- 
formation, as 



(4) 



An overlap between different mo- and mi-quasiparticle 
states 



(o) 
fei,- 



(i) 



(* (0) l^---/3S /3«t...^t|$( 1 ) )) (5 ) 



is decomposed to a linear combination of the following 
matrix elements with 2m = uiq + mi, 



Mi = M {l 

<*<°>|Cfc... 



Jk f 



,t 



4 

J2r 



|$(D) 



(6) 



where 1 < jk < M, 1 < j' k < M and < k < 2m. Indices 
jk and j' k are relabelled with i's as ii < 12 < • • • < 
%2m- Hereafter, we denote this set of indices {ik',k — 
1, 2, 3, • • • , 2m} as /. / is a subset of [2M], that is, I = 
{ii,*a,--- ,»am} C {1,2,3,- •• ,2M} = [2M]. Here we 
use a compact notation for a group of indices by following 
the Ref.[15|. This notation gives us a good view for a 
formula with complicated indices. 

It is convenient to introduce a new expression di for 



di = 



(i = l,-..,M) 

(i = 1 + M, ■ ■ ■ , 2M) 



(7) 



where i = 2M + 1 — i. By this definition, an ordering 
concerning the d operators in Eq.© can be systemati- 
cally handled. As a result, we can deal with products 
consisting of the arbitrary number of c and in a uni- 
fied way. It should be noted that in the ordering of the 
creation-annihilation operators in Eq.([5]), all the annihi- 
lation operators c are placed in the left of the creation 
operators c'. This arrangement needs to be considered 
for the convenience of the mapping to Grassmann num- 
bers. 



III. MAPPING BY FERMION COHERENT 
STATE 

By introducing the Fermion coherent state[l[ 

|O = e-£.€"4| >, (8) 

we map the creation and annihilation operators to Grass- 
mann numbers £* and £ respectively, which follow the 
anticommutation rules 



Sid + 66 



(9) 
(10) 
(11) 



By definition, the coherent state is introduced as an 
eigenstate of the annihilation operator, 



ciie> = &io. 



(12) 



where & needs to be a Grassmann number due to the an- 
ticommutation nature of the annihilation operators. The 
conjugate relation is 



(13) 
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where £* is a conjugate variable to £j. An expectation 
value of the d operator in Eq. Q is mapped to a Grass- 
mann number as, 



(Z\dAt) = 6, 



(14) 



where £j corresponds to £* for i = 1, ■ • • , M and £g for 
i = 1 + M, ■ • • , 2M. A vector d consisting of operators c 
and c' as 

d = ( C I> C L--- ,c\ l7 c M ,c M -ir ■ ■ »ci) (15) 
can be mapped to a Grassmann vector £ as 

f = (££,•••,&,&,,&,-!,••• (i6) 

Next, we map the general matrix element Mi in Eq.© 
to Grassmann numbers by inserting the completeness 



(17) 



between Cj h and cX, in Eq.©. By noting that (0|£) = 1, 
we can obtain an expression of the matrix element Mi 
in Eq.® in terms of the following Grassmann integral 



Mi= I P£exp( -eX|)& 2 „ 



■k, (18) 



where V£ = d£ 2M • • • d& = (-1) M Hereafter 
M is assumed to be even. The matrix element Mi can 
be shown by the Grassmann integral in terms of only 
£. This fact directly leads to a Pfaffian form as will be 
discussed in a next section. 

In the zero quasiparticle case, the corresponding ma- 
trix element is just a norm-overlap kernel of the original 
HFB state Eq.©. Robledo @ was the first to give a 
proof that the norm-overlap is represented in terms of 
the Pfaffian, that is, 



($(o)|^(i)) 



Pf(X), 



(19) 



where a convention in Ref . (Tol j is employed for Grassmann 
vectors in this expression. This convention is different 
from Ref. [2] in the ordering of Grassmann numbers in 
Eq. (I16p . With this convention, a matrix X is given as 



Z« -A 
A -AZ(°)*A 



where A is defined as 
A 



(20) 



(21) 



Thanks to the Pfaffian in the formula, an ambiguity in 
the sign assignment is removed, which was the primary 
concern in the Onishi formula having an expression of a 
square root of a determinant. Note that, by Eq. ([2"5j) in 



the next section, we can immediately obtain the Pfaffian 
form for the norm-overlap kernel Ea. (ll9D . 

As for matrix elements with respect to one- 
quasiparticle HFB states, we have derived a formula 
in Ref. 10]. With this formula, quantities such 
as ($(°)| Cl cj|$«}, (fcWlcJctlsW), and ($W\ac j \& 1 '>), 
can be evaluated. In our proof, a variable transforma- 
tion fj — is introduced. A factorization of the anti- 
symmetric matrix X is carried out with a non-singular 
matrix R, as 



X = i?*XR, 



(22) 



so as to derive the formula. The matrix J in the above 
expression is given as 



J 



/ 
-I 



(23) 



where / is the unit matrix with the M dimension. The 
resultant formula is given in terms of 1 



(&°')\d j di\& 1 ')) = ($ (0) |$ (1) )(X 



as 



(-1)1 



= P 



(24) 



The formula shown in Ref. [Toj ] contains extra matrix ele- 
ments consisting of the U and V factors, but in the above 
expression we show only the essential part so as to be 
consistent with the generalization to multi-quasiparticle 
states, to be discussed below. 

In order to derive a new formula for general matrix 
elements, we consider a generating function with Grass- 
mann numbers in the subsequent sections. 



IV. GRASSMANN INTEGRALS 

In this section, we review the Pfaffian type and deter- 
minant type of Grassmann integrals. By cosidering these 
two types by comparison, it is possible to understand 
how natural the Pfaffian is in handling the HFB theory 
in a view of mathematical structure. Based on the Refs. 
0, H^h w e emphasize that the Pfaffian is mathemati- 
cally more fundamental, and that it plays an important 
role in the calculation of the matrix elements in the HFB 
theory. Here, we use #'s and z's instead of £'s, for the 
sake of generality. 

Firstly, the Pfaffian-type Grassmann integral is given 

by 



/• 



d9 2n ■■■dQ x exp -FAB = Pf(A), (25) 



where 6\ , 6%, ■ ■ ■ , 02n are Grassmann variables and 8 f = 
(0i)02) - "" ,#2n) is a Grassmann vector. Matrix A is a 
2n x 2n skew-symmetric matrix. It can be easily proved 
by expanding the exponent, as shown in Refs.[16l.ll7j. for 
instance. 
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by 



The determinant type of Grassmann integral is given 



for all k and I G /. From a viewpoint of Wick's 
theorem, this relation implies that the matrix element 
' ' ' z i m %i • • • Z ?) can be rewritten in terms of the 



Yl dz* a dz a exp (-z^Bz) = det(B) (26) contractions (ziZ*)'s because of (B — (z;z*). This 



where z±, Z%, • • • ,z n are Grassmann variables and z* 
is a conjugate Grassmann number of Zj. The z 4 = 
(zi,Z2,--- ,z n ) is a Grassmann vector and and z* is a 
conjugate vector of z. The B is an n x n matrix. Its 
proof is given, for instance, in Rcf.[l]. 

Depending on how to define a Grassmann vector and 
how to set a form of the integrand, the outcome of the 
integral changes: it takes a form of either a determinant 
or the Pfaffian. It is worth noting, however, that the de- 
tcrminat type of Grassmann integral is included as a spe- 
cial case of the Pfaffian- type Grassmann integral [TtJ ■ To 
demonstrate this fact, let us consider a skew-symmetric 
matrix A in the following form. 



A = 





(BAY 



-BA 




(27) 



where B is the matrix appearing in Ea. (|26[) . In the 
Grassmann vector in Ea. (j25|) . we choose 9% = z* for 
i = 1, • • • , n and 9i — Z2n+i-i for i = n + 1, • • • , 2n, 
namely, a bipartite representation 6 t = (z*,Az) is em- 



ployed here, 
of 



and 



Then, Ea.([23|l reduces to Eg. (121)1) because 



Pf(A) = (-l) n det(B), (28) 



d6 2n ---d6 1 = (-l) n Y[dz* a dz a . (29) 

o*=l 

From this result, it is possible to regard that the Pfaffian- 
type Grassmann integral is more fundamental. We call 
this property the determinant-Pfaffian correspondence. 

The Grassmann integrals and Wick's theorem are 
closely related to each other. It is known that a 
technique of generating functions is useful for finding 
a relation between them By making use of the 

determinant-type Grassmann integral Eq. (l26[) . matrix el- 
ement {zi t • ■ ■ Zi m z* m ■ ■ ■ z* } can be expressed in a com- 
pact way as 

{ z i! ' ' ' z i m z i m ' ' ' z i 1 ) 

I IL dz* a dz a exp {\z* l Bz) z h ■ ■ ■ z i m z * m ■■■z* 1 - 



I Ua d z adz a exp {\z* l Bz) 



det{{B- 1 )!), 



(30) 



with m < n. This relation can be obtained with a gen- 
erating function proposed in Ref.[l|. As introduced in 
Sec. II, a compact notation for a group of indices fl5| is 
employed here again. An mxm sub-matrix (B~ 1 )j with 
I = {ii,i2, • ■ • , im}, is defined as 



result is general and needs no assumption of specific sys- 
tems for the derivation, which is presented in Eqs.(2.84)- 
(2-87) of Ref.d 

through a generating function. 
On the contrary, for the evaluation of the HFB matrix 
elements Ea. ([T8"|l . the Pfaffian version of Ea. ([30l is useful. 
In a similar manner to the determinant case, we expect 
that (£i 2m • • ■ ) should be decomposed into combina- 
tions of the contractions (£i£j)'s. The existence of this 
relation is likely due to the determinant-Pfaffian corre- 
spondence. Its derivation can, in fact, be easily carried 
out in the same procedure with a generating function, 
which will be discussed in the next section in detail. We 
learned that this Pfaffian relation is briefly mentioned in 
Ref . [16| , and it is summarized in appendix A of a very re- 
cent paper [l8| from a mathematical point of view. How- 
ever, in the context of a physical application of the rela- 
tion, we are the first to present the relation in detail how 
to evaluate the matrix elements. In addition, a compar- 
ison of the result to the conventional expansion by the 
usual Wick theorem is demonstrated for the first time in 
this paper. 



V. GENERATING FUNCTION 

To evaluate general matrix elements for multi- 
quasiparticle states, we introduce a generating function 
defined by 

G(J) ee ^— - — y 2 - s 2 I. 2S ; 32 

where a "source" J is defined as a Grassmann vector 
J ee (J*, Jj, • • • , Jm, Jm, Jm-i, ■ ■ ■ , J\)- (33) 
Introducing a new variable £' by shifting £ as 

the following relation is obtained 



(34) 



222 2 



-JtX^J. 
2 



The generating function thus becomes 



G(J) = cxp ( |J*X- X J 



(35) 



(36) 



{{B~ l )i)k,i = [B~ l )i h 



(31) 



A differentiation of the original form of the generat- 
ing function Eg . ((32)) with respect to the source vector J, 
gives rise to the matrix element Aij in Eg. (fig)) in the 



5 



form of 



VI. GENERALIZED WICK'S THEOREM 



5< 



2m) 



G(J) 



/2?£exp (K'Xe + PC - |l* J) & 



/Pfexp(^X£) 
/^exp(i^4 m ---4 



j=o 



J^exp(ie 



(37) 



On the other hand, by differentiating the form of 
Eg. p6p . we can obtain an alternative expression for 



S (2* 



■5./, 



-G(J) 



j=o 



For the aim to express in a compact 



manner, we introduce a 2m-dimensional sub-vector Jj 
and a 2m x 2m sub-matrix Aj, with J = {ij, i 2 , • • • , i 2m }- 
They are defined as 



(J/)fe 
(^i)fc,l 



= -A, 
= j4. 



(38) 



for all A; and I € I. As the Grassmann numbers in a sub- 
group [2M] — I vanishes, the derivative of the generating 
function is calculated to be 



SJi. 



G{J) 



j=o 



6 J. 



A( 2m ) /I 

-^-exp(-J/(X-^ 



Pf((x- 1 ) / ). 



(39) 



In the above calculation, we use a useful relation for any 
skew-symmetric matrix A with 2m x 2m dimension and 
2m Grassmann numbers z's, which is 



Sz 2m ■■■8zi 



exp (A) 
1 . ™ 



Sz 2m ■■■5zi ml 

g(2m) 

PKA) Zl Z 2 ■ 



Sz 2m ■ 

= PUl 

1 v~v2m 
2 



■ Szi 



■ z 2r , 



(40) 



where A = ^ VJ* "j—i z\Ai^ Z j. Eo. (j52l) in Appendix is 
also used for the above derivation. 

Equating Eqs. (|3"Tl) and (I3T)1) . a new formula is obtained 
for matrix elements between different multi-quasiparticle 
states: 



Mi = (*W|*( 1 ))Pf((X- 1 ) J ), 
= Pf(X)Pf((X- 1 ) / ). 



(41) 



This is a generalization of Eg . (|2~4l because Pf(A) = a 
as in Eg. (|53[) for a 2 x 2 anti-symmetric matrix A = 
a 



-a 



It should be noted here that we can also 



derive Eg. ([3"Ul) from Eg. (l4"TT) by the dctcrminant-Pfaffian 
correspondence, which is to be presented in Section IV. 

As we will discuss below, our formula turns out to be an 
extension of the traditional generalized Wick theorem^. 



The standard approach to calculate matrix elements 
of many-body operators is to use generalized Wick's the- 
orem [3j, in which a chain of operators is expanded in 
terms of normal orders and contractions. 

Contractions in the HFB theory correspond to the 
transition densities (p^ 01 ' and p' 01 -*), and the transition 
pairing tensors (k^ 01 ) and k^ 10 )*). These quantities are 
mathematically expressed as 1 



(*W\c]ci\$W) 

($(0)|$(1)} ' 

($(°)| Ci 4|$W) 

($(0)|$(1)) ' 
Ki >3 ($(0)|$(1)) ' 

,(io)* _ ^ (0) l4c]|$W) 



(01) 

Pij 



($(0)|$(1)) 



(42) 



where 1 < i,j < M. In accordance with generalized 
Wick's theorem, any matrix elements in Eq.([B]) can be 
expressed by these contractions in a combinatorial way. 

In the present method, these basic contractions are 
connected to X -1 directly, for instance, 



(10). 



-(*(°>|dtdt|* (1) ) 

($(0)|$(1)) 



(43) 



where Eg. (|2"4")l is used. In general, a relation between 
X _1 ,/9, and k is given as 



X" 



_ K (io)* ^(01)1 
_p(01) K (oi) 

_ K (10)* (l_ p (01))t 



,(oi) _ i 



H 



(01) 



(44) 



The transition densities and transition pairing tensors 
can be also expressed in terms of the skew-symmetric 
matrices in the Thouless form Eq.([T]), as shown in 
Appendix E of Ref.[3j, which needs somewhat lengthy 
calculations for the proof. In the present method, how- 
ever, X -1 is directly calculated from Eo. (f2T))) . which re- 
sults in 



1 
A 



1 
A 



(45) 



I (zwz*w - zW(i -z<°yzwy 



1 In Ref. p( 01 ' is employed as the transition density because 
p( 01 ) = 1 — p' 01 ', but we use p( 01 ) in the present work. 
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This result is consistent with the expressions (E.54) in 
Ref.Q. 

Let us now confirm that our formula to calculate many- 
body matrix elements is consistent with the result ob- 
tained with generalized Wick's theorem. As an exam- 
ple, let us consider a case of two-quasiparticle states 
|$^ 2 ) = /3l 1 Pl 2 \^ p) ). In the evaluation of the rele- 
vant norm-overlap, several types of many-body matrix el- 
ements need to be calculated, say (<&(°)|c; 4 c; 3 c| q 
This quantity can be calculated by means of our formula 
Eq.(gIJ as, 



($ (0) |Q 4 Q3«|$(i)) 



(46) 



($(o)|$(i)) 

= Pfttx- 1 )^^) 
+(x- 1 ) (il4 (x- 1 ) i2j3 . 

- _ K (10) *K (01) - n m n m + n {01) n {01) 

where < h,h,l3,h < M. In the above equations, the 
second line is obtained with the formula Eq. (|4"Tj) , and the 
third line is an explicit expansion of the Pfaffian due to 
the definition Eq. ([54")) . The last line is a substitution with 
the result in Ea. (|44|) . The last line is consistent with the 
result derived from traditional Wick's theorem. In other 
words, the present Pfaffian formula is an alternative rep- 
resentation of Wick's theorem. For details of generalized 
Wick's theorem, see Ref.[4J for instance. 

Recently, Perez-Martin and Robledo also reformu- 
lated Wick's theorem based on the statistical mechanics 
(Gaudin's theorem) to give a closed expression in Ref. 
11]. Also, Bcrtsch and Robledo presented another closed 
form with the Pfaffian [T^]. In comparison with these for- 
mulae, the present formula has a fairly compact and sim- 
ple expression. The main reason for this compactness 
and simplicity comes from an introduction of X" 1 , by 
which the transition density and transition pairing ten- 
sor can be unified. As a consequence, the expression of 
Aii becomes free from combinatorial complexity in ap- 
plying generalized Wick's theorem. 



VII. BALIAN-BREZIN CONJECTURE 

In nuclear structure physics, generalized Wick's the- 
orem is quite important especially for a description of 
multi-quasiparticle states based on an HFB vacuum. 
Here we revisit the contents of the theorem historically. 

In Ref. [13] , Lowdin first noticed that an identity for de- 
terminants, called the Lewis Carroll (Dodgson) formula 
or Desnanot-Jacobi adjoint matrix theorem [19j. can be 
applied to the density matrix for Slater determinants [l4j] . 



In Ref. 14], Balian and Brezin conjectured that a Pfaffian 



counterpart of Sylvester's identity[2l[ for Slater deter- 
minants would be useful for the HFB theory. However, 
there has been no progress in this conjecture since its pro- 



posal in 1969 as shown in the first footnote of Ref. [131]. 
In this section, we discuss this conjecture. 

To approach the Balian-Brezin conjecture, we consider 
the Ref. [131 ] where they consider an another way to carry 
out the Grassmann integrations in Eq. (fT51) . The Grass- 
mann variables {£J , Q , • • • , £* M , £ M , £ M - 1 , • • • , £i } are di- 
vided into two groups. The first group of the Grass- 
mann variables consists of I = • ■ • , i2m}, which cor- 
responds to the indices appearing in the operators di in 
Eq.©. The other group consists of the rest of the indices 
in the total model space, which is denoted as / = [2M] — I 
meaning the complimentary group of /. 

Due to a reordering of the Grassmann variables, the 
integral measure is rewritten as 



IJd£de« = Z>&2>&(-) 



\i\ 



(47) 



and 



where \I\_ = J2k=i % Here = ' ' ' d^. 

V£j = d!; j2M _ 2m ■ ■ ■ d£ jl2 [ji < j 2 < ■ ■ ■ < 32M-2m,jk G 

/). Integration over X>£j is easily carried out and gives 
rise to unity. The matrix element for A4j thus becomes 



Mt = (-)l'l 



(48) 



Here, we use the notation Aj for a 2(M — m) x 2(M — m) 
sub-matrix of a matrix A, which is produced by remov- 



ing the matrix elements Ai j for i,j 



,l2r, 



fror 



the original matrix A. Note that the resultant formula 
Eq. (|48p is nothing to do with generalized Wick's theorem, 
unlike the formula we have derived in Eq. (|41[) . 

Nonetheless, it is possible to relate our result with the 
above result, through an identity for the Pfaffian. By 
comparing Eq. (|4"Tj) to Eq. (l4"51) . there should be a relation 



(49) 



Pf(X)Pf((X" 1 ) / ) = (-) |7| Pf(X r ). 



We found that this relation holds for any skew-symmetric 
matrix X, and that this Pfaffian identity has been re- 
cently discovered by mathematicians [l^]. Their proof is, 
however, given in an elementary way without relying on 
any Grassmann integrations [la . |22j. This identity is 
called the Pfaffian version of the Lewis Carroll (Dodg- 
son) formula [15| , whose original version holds for deter- 
minant. 

Through the comparison between our result and the 
one obtained by Avez and Bender [13| . we found that 
the both Pfaffian formula are firmly connected by the 
Pfaffian identity, which is just what Balian and Brezin 
conjectured for the HFB caseflj]. 



VIII. SUMMARY 

In the present paper, we presented a compact and 
closed formula to evaluate matrix elements for multi- 
quasiparticle HFB states, and a relation to generalized 
Wick's theorem is discussed. 
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To calculate multi-quasiparticle matrix elements, we 
started our derivation from a fact that the Thouless 
form of the HFB wave function can be rewritten by the 
Fermion coherent state with Grassmann numbers, by fol- 
lowing Robledo's approach [2(]. By a generating function 
with Grassmann numbers, we could derive a new Pfaf- 
han formula to calculate matrix elements between multi- 
quasiparticle HFB states. This formula is an extension 
of the norm-overlap kernels of one-quasiparticle states 
which we presented in Ref. [10]. In traditional general- 
ized Wick's theorem, the transitional density and pairing 
tensor are utilized in a combinatorial way, while, in this 
new formula, these basic contractions are utilized with no 
distinction. In consequence, the formula becomes simple 
and compact and is free from combinatorial complexity 

Besides, by comparing our Pfaffian formula to the 
one obtained by Avez and Bender [l3j], we found that 
the both Pfafnan formulae are closely related to the 
Pfafnan version of the Lewis Carroll formula in pure 
mathematics [Pol] . In Ref . [blj . Balian and Brezin con- 
jectured a usefulness of extended Sylvester's identity by 
means of the Pfafnan [2l| . The present paper gives an 
affirmative answer to their conjecture. 

Appendix 

For a skew-symmetric matrix A with dimension 2n x 
2n, whose matrix elements are a%j , the Pfafnan is defined 



as 

1 ™ 
Fi{A) =2^ ^ Bgn(«7)JJa a(aj _i )<r(2i) (50) 
' ctsS 2 „ i=l 

where a is a permutation of {1,2, 3,- •■ , 2n}, sgn(<r) is 
its sign, and S^n represents a symmetry group. 

For Grassmann numbers Zi, z%, ■ ■ ■ , Z2„, if we define Z 

as 

Z = y^ajjZjZj, (51) 

then Z n is shown by the Pfaffian as 

-Z n = Vi{A)z lZ2 ---z 2n . (52) 

This is a definition of the Pfaffian in the exterior algebra. 

According to Ea. (|50p . for a n x n (n = odd) matrix, 
Pi(A) = 0. 

For a 2 x 2 matrix, 

Pf(A)=oi 2 . (53) 

For a 4 x 4 matrix, 

Pf(A) = ai 2 a 34 - a u a 2 4 + aua 2 z. (54) 
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